Microbial Succession in the Gut: Directional Trends of Taxonomic and Functional Change in a Birth Cohort of Spanish Infants

In spite of its major impact on life-long health, the process of microbial succession in the gut of infants remains poorly understood. Here, we analyze the patterns of taxonomic and functional change in the gut microbiota during the first year of life for a birth cohort of 13 infants. We detect that individual instances of gut colonization vary in the temporal dynamics of microbiota richness, diversity, and composition at both functional and taxonomic levels. Nevertheless, trends discernible in a majority of infants indicate that gut colonization occurs in two distinct phases of succession, separated by the introduction of solid foods to the diet. This change in resource availability causes a sharp decrease in the taxonomic richness of the microbiota due to the loss of rare taxa (p = 2.06e-9), although the number of core genera shared by all infants increases substantially. Moreover, although the gut microbial succession is not strictly deterministic, we detect an overarching directionality of change through time towards the taxonomic and functional composition of the maternal microbiota. Succession is however not complete by the one year mark, as significant differences remain between one-year-olds and their mothers in terms of taxonomic (p = 0.009) and functional (p = 0.004) microbiota composition, and in taxonomic richness (p = 2.76e-37) and diversity (p = 0.016). Our results also indicate that the taxonomic composition of the microbiota shapes its functional capacities. Therefore, the observed inter-individual variability in taxonomic composition during succession is not fully compensated by functional equivalence among bacterial genera and may have important physiological consequences. Finally, network analyses suggest that positive interactions among core genera during community assembly contribute to ensure their permanence within the gut, and highlight an expansion of complexity in the interactions network as the core of taxa shared by all infants grows following the introduction of solid foods

Analysis of among-site variation in substitution patterns

Substitution patterns among nucleotides are often assumed to be constant in phylogenetic analyses. Although variation in the average rate of substitution among sites is commonly accounted for, variation in the relative rates of specific types of substitution is not. Here, we review details of methodologies used for detecting and analyzing differences in substitution processes among predefined groups of sites. We describe how such analyses can be performed using existing phylogenetic tools, and discuss how new phylogenetic analysis tools we have recently developed can be used to provide more detailed and sensitive analyses, including study of the evolution of mutation and substitution processes. As an example we consider the mitochondrial genome, for which two types of transition deaminations (C⇒T and A⇒G) are strongly affected by single-strandedness during replication, resulting in a strand asymmetric mutation process. Since time spent single-stranded varies along the mitochondrial genome, their differential mutational response results in very different substitution patterns in different regions of the genome

Evolutionary Analysis of Mitogenomes from Parasitic and Free-Living Flatworms

Copyright: © 2015 Solà et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. The attached file is the published version of the article

Comprehensive analysis of the base composition around the transcription start site in Metazoa

BACKGROUND: The transcription start site of a metazoan gene remains poorly understood, mostly because there is no clear signal present in all genes. Now that several sequenced metazoan genomes have been annotated, we have been able to compare the base composition around the transcription start site for all annotated genes across multiple genomes. RESULTS: The most prominent feature in the base compositions is a significant local variation in G+C content over a large region around the transcription start site. The change is present in all animal phyla but the extent of variation is different between distinct classes of vertebrates, and the shape of the variation is completely different between vertebrates and arthropods. Furthermore, the height of the variation correlates with CpG frequencies in vertebrates but not in invertebrates and it also correlates with gene expression, especially in mammals. We also detect GC and AT skews in all clades (where %G is not equal to %C or %A is not equal to %T respectively) but these occur in a more confined region around the transcription start site and in the coding region. CONCLUSIONS: The dramatic changes in nucleotide composition in humans are a consequence of CpG nucleotide frequencies and of gene expression, the changes in Fugu could point to primordial CpG islands, and the changes in the fly are of a totally different kind and unrelated to dinucleotide frequencies

Shifting patterns of natural variation in the nuclear genome of caenorhabditis elegans

<p>Abstract</p> <p>Background</p> <p>Genome wide analysis of variation within a species can reveal the evolution of fundamental biological processes such as mutation, recombination, and natural selection. We compare genome wide sequence differences between two independent isolates of the nematode <it>Caenorhabditis elegans </it>(CB4856 and CB4858) and the reference genome (N2).</p> <p>Results</p> <p>The base substitution pattern when comparing N2 against CB4858 reveals a transition over transversion bias (1.32:1) that is not present in CB4856. In CB4856, there is a significant bias in the direction of base substitution. The frequency of A or T bases in N2 that are G or C bases in CB4856 outnumber the opposite frequencies for transitions as well as transversions. These differences were not observed in the N2/CB4858 comparison. Similarly, we observed a strong bias for deletions over insertions in CB4856 (1.44: 1) that is not present in CB4858. In both CB4856 and CB4858, there is a significant correlation between SNP rate and recombination rate on the autosomes but not on the X chromosome. Furthermore, we identified numerous significant hotspots of variation in the CB4856-N2 comparison.</p> <p>In both CB4856 and CB4858, based on a measure of the strength of selection (k_a/k_s), all the chromosomes are under negative selection and in CB4856, there is no difference in the strength of natural selection in either the autosomes versus X or between any of the chromosomes. By contrast, in CB4858, k_a/k_svalues are smaller in the autosomes than in the X chromosome. In addition, in CB4858, k_a/k_svalues differ between chromosomes.</p> <p>Conclusions</p> <p>The clear bias of deletions over insertions in CB4856 suggests that either the CB4856 genome is becoming smaller or the N2 genome is getting larger. We hypothesize the hotspots found represent alleles that are shared between CB4856 and CB4858 but not N2. Because the k_a/k_sratio in the X chromosome is higher than the autosomes on average in CB4858, purifying selection is reduced on the X chromosome.</p

Evolution and Dynamics of Regulatory Architectures Controlling Polymyxin B Resistance in Enteric Bacteria

Complex genetic networks consist of structural modules that determine the levels and timing of a cellular response. While the functional properties of the regulatory architectures that make up these modules have been extensively studied, the evolutionary history of regulatory architectures has remained largely unexplored. Here, we investigate the transition between direct and indirect regulatory pathways governing inducible resistance to the antibiotic polymyxin B in enteric bacteria. We identify a novel regulatory architecture—designated feedforward connector loop—that relies on a regulatory protein that connects signal transduction systems post-translationally, allowing one system to respond to a signal activating another system. The feedforward connector loop is characterized by rapid activation, slow deactivation, and elevated mRNA expression levels in comparison with the direct regulation circuit. Our results suggest that, both functionally and evolutionarily, the feedforward connector loop is the transitional stage between direct transcriptional control and indirect regulation

Atypical AT Skew in Firmicute Genomes Results from Selection and Not from Mutation

The second parity rule states that, if there is no bias in mutation or selection, then within each strand of DNA complementary bases are present at approximately equal frequencies. In bacteria, however, there is commonly an excess of G (over C) and, to a lesser extent, T (over A) in the replicatory leading strand. The low G+C Firmicutes, such as Staphylococcus aureus, are unusual in displaying an excess of A over T on the leading strand. As mutation has been established as a major force in the generation of such skews across various bacterial taxa, this anomaly has been assumed to reflect unusual mutation biases in Firmicute genomes. Here we show that this is not the case and that mutation bias does not explain the atypical AT skew seen in S. aureus. First, recently arisen intergenic SNPs predict the classical replication-derived equilibrium enrichment of T relative to A, contrary to what is observed. Second, sites predicted to be under weak purifying selection display only weak AT skew. Third, AT skew is primarily associated with largely non-synonymous first and second codon sites and is seen with respect to their sense direction, not which replicating strand they lie on. The atypical AT skew we show to be a consequence of the strong bias for genes to be co-oriented with the replicating fork, coupled with the selective avoidance of both stop codons and costly amino acids, which tend to have T-rich codons. That intergenic sequence has more A than T, while at mutational equilibrium a preponderance of T is expected, points to a possible further unresolved selective source of skew

Bacterial Niche-Specific Genome Expansion Is Coupled with Highly Frequent Gene Disruptions in Deep-Sea Sediments

The complexity and dynamics of microbial metagenomes may be evaluated by genome size, gene duplication and the disruption rate between lineages. In this study, we pyrosequenced the metagenomes of microbes obtained from the brine and sediment of a deep-sea brine pool in the Red Sea to explore the possible genomic adaptations of the microbes in response to environmental changes. The microbes from the brine and sediments (both surface and deep layers) of the Atlantis II Deep brine pool had similar communities whereas the effective genome size varied from 7.4 Mb in the brine to more than 9 Mb in the sediment. This genome expansion in the sediment samples was due to gene duplication as evidenced by enrichment of the homologs. The duplicated genes were highly disrupted, on average by 47.6% and 70% for the surface and deep layers of the Atlantis II Deep sediment samples, respectively. The disruptive effects appeared to be mainly due to point mutations and frameshifts. In contrast, the homologs from the Atlantis II Deep brine sample were highly conserved and they maintained relatively small copy numbers. Likely, the adaptation of the microbes in the sediments was coupled with pseudogenizations and possibly functional diversifications of the paralogs in the expanded genomes. The maintenance of the pseudogenes in the large genomes is discussed

The enigmatic mitochondrial genome of Rhabdopleura compacta (Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

<p>Abstract</p> <p>Background</p> <p>The Hemichordata comprises solitary-living Enteropneusta and colonial-living Pterobranchia, sharing morphological features with both Chordata and Echinodermata. Despite their key role for understanding deuterostome evolution, hemichordate phylogeny is controversial and only few molecular data are available for phylogenetic analysis. Furthermore, mitochondrial sequences are completely lacking for pterobranchs. Therefore, we determined and analyzed the complete mitochondrial genome of the pterobranch <it>Rhabdopleura compacta </it>to elucidate deuterostome evolution. Thereby, we also gained important insights in mitochondrial tRNA evolution.</p> <p>Results</p> <p>The mitochondrial DNA of <it>Rhabdopleura compacta </it>corresponds in size and gene content to typical mitochondrial genomes of metazoans, but shows the strongest known strand-specific mutational bias in the nucleotide composition among deuterostomes with a very GT-rich main-coding strand. The order of the protein-coding genes in <it>R. compacta </it>is similar to that of the deuterostome ground pattern. However, the protein-coding genes have been highly affected by a strand-specific mutational pressure showing unusual codon frequency and amino acid composition. This composition caused extremely long branches in phylogenetic analyses. The unusual codon frequency points to a selection pressure on the tRNA translation system to codon-anticodon sequences of highest versatility instead of showing adaptations in anticodon sequences to the most frequent codons. Furthermore, an assignment of the codon AGG to Lysine has been detected in the mitochondrial genome of <it>R. compacta</it>, which is otherwise observed only in the mitogenomes of some arthropods. The genomes of these arthropods do not have such a strong strand-specific bias as found in <it>R. compacta </it>but possess an identical mutation in the anticodon sequence of the tRNA_Lys.</p> <p>Conclusion</p> <p>A strong reversed asymmetrical mutational constraint in the mitochondrial genome of <it>Rhabdopleura compacta </it>may have arisen by an inversion of the replication direction and adaptation to this bias in the protein sequences leading to an enigmatic mitochondrial genome. Although, phylogenetic analyses of protein coding sequences are hampered, features of the tRNA system of <it>R. compacta </it>support the monophyly of Ambulacraria. The identical reassignment of AGG to Lysine in two distinct groups may have occurred by convergent evolution in the anticodon sequence of the tRNA_Lys.</p

The Evolution of Mammalian Gene Families

Gene families are groups of homologous genes that are likely to have highly similar functions. Differences in family size due to lineage-specific gene duplication and gene loss may provide clues to the evolutionary forces that have shaped mammalian genomes. Here we analyze the gene families contained within the whole genomes of human, chimpanzee, mouse, rat, and dog. In total we find that more than half of the 9,990 families present in the mammalian common ancestor have either expanded or contracted along at least one lineage. Additionally, we find that a large number of families are completely lost from one or more mammalian genomes, and a similar number of gene families have arisen subsequent to the mammalian common ancestor. Along the lineage leading to modern humans we infer the gain of 689 genes and the loss of 86 genes since the split from chimpanzees, including changes likely driven by adaptive natural selection. Our results imply that humans and chimpanzees differ by at least 6% (1,418 of 22,000 genes) in their complement of genes, which stands in stark contrast to the oft-cited 1.5% difference between orthologous nucleotide sequences. This genomic “revolving door” of gene gain and loss represents a large number of genetic differences separating humans from our closest relatives